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Data from the Kepler satellite (Q0-Q1 1) are used to study HAT-P-7. The satellite's data are extremely valuable for astero- 
seismic studies of stars and for observing planetary transits; in this work we do both. An a steroseismic study of the host star 
improves the accuracy of the stellar parameters derived by Christensen-Dalsgaard et al. (2010), who followed largely the 
same procedure but based the analysis on only one month of Kepler data. The stellar information is combined with transit 
observations, phase variations and occultations to derive planetary parameters. In particular, we confirm the presence of 
ellipsoidal variations as discovered by Welsh et al. (2010 ), but revise their magnitude, and we revise the occultation depth 
( [Boruck i et al. 2009 ), which leads to different planetary temperature estimates. All other stellar and planetary parameters 
are now more accurately determined. 
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1 Introduction 

The Kepler mission was launched on March 6, 2009. It was 
designed explicitly to be able to detect earth- size planets, 
including those in the habitable zone. The detection method 
for Kepler to achieve this goal is measuring planetary tran- 
sits to a high level of photometric precision. To improve 
the chances of detection, a large number of stars have to 
be monitored for a long period of time, preferably continu- 
ously. It was quickly realised that Kepler would also provide 
extremely valuable data to study stellar pulsations. Data for 
this purpose are obtained from a subset of the about 500 
stars that are sampled every 58.8 seconds (short cadence), 
out of a much larger set of about 150 000 stars observed at 
a rate of 29.4 minutes ( [Borucki et al.|2008| ). 

The satellite orbits in an earth-trailing heliocentric orbit 
(ETHO), with a 55° Sun avoidance angle, and continuously 
observes a patch of the sky centred on the Cy gnus-Lyra re- 
gion. The data are organised into quarters; after each quarter 
the photometer is rolled 90° to keep the solar arrays pointed 



et al.| p010| ), based on early Kepler data (Q0-Q1). After 



at the Sun. This causes a 42 hour data gap ( Haas et al. 2010} . 
This work uses Q0-Q1 1 data of HAT-P-7, which has been 
continuously observed in Kepler's short cadence mode. 

HAT-P-7 (KIC 10666592, visual magnitude 10.46, spec- 
tral type F8) has been identified as a planet host by ground- 
based observations prior to the start of the Kepler mission 
( Pal et al.]|2008 ). It is one of the most extensively studied 
star-planet systems known to date. A study of the stellar 
pulsations has been carried out by Christensen-Dalsgaard] 
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the initial discovery of its exoplanet HAT-P-7b ( |Pal et al. 
2008), phase variations and planetary occultations were dis 
covered in Kepler's Q0 data ( [Borucki et al.|20 09), and con- 
sequently the planetary temperature and atmosphere were 
studied ( [Christiansen et al.|2010||Spiegel & Burrows|2010| ). 
Additional ellipsoidal variations were later discovered in the 
light curve ( We lsh et~aL]|2010| ) and the system is currently 
the only planetary case where this effect has been found. 
HAT-P-7 has a surprising inclination angle: a study of the 
Rossiter-McLaughlin effect in radial velocity data has found 
it to be retrograde or polar ( Winn et al. 2009| ). These authors 
also suggest the presence of a third body. A study of tran- 
sit timings based on EPOXI data found no evidence of an 
additional planet ( [Ballard et al.|2011| ). 

Here, we make a preliminary anaysis of the system 
based on almost three years of Kepler data. We derive stellar 
pulsation frequencies and perform an asteroseismic mod- 
elling to derive the properties of the host star. The plane- 
tary system is then studied making use of these newly de- 
termined stellar properties. 

2 Stellar properties 

Stellar pulsations cause temperature variations, which lead 
to flux variations. The frequencies are observed from the 
power spectrum. Different models are then compared with 
the observed frequencies and the best fit leads to the stel- 
lar properties. The procedure followed here resembles the 
asteroseismic analysis of HAT-P-7 based on only the first 
month of Kepler data ( Christensen-Dalsgaard et al.|2010 ). 
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Fig. 1 Top figure: oversampled (red) and smoothed power 
spectrum (black), showing regular spacings governed by the 
asymptotic equation (eq. [T}. Bottom figure (zoom): mode 
identification for three observed frequencies. An algorithm 
determines the maxima in the smoothed spectrum, the rele- 
vant maxima (vertical lines) are selected manually making 
use of the asymptotic equation. Sometimes, an average of 
two maxima is used (e.g. I = 2 mode in this figure). 



2.1 Observations and data analysis 

The raw Kepler data are detrended and normalised, and 
transit-like features are removed using a specifically de- 
signed median filter. The filtered photometry is then used 
to calculate a power spectrum. A spectrum with a sampling 
rate of 0.005 /iHz is used, which means an oversampling of 
about three times. A median filter with a period of 1 /iHz 
is used to smooth the power spectrum, to remove the fine 
structure of the spectrum, caused by the finite lifetime of 
the modes. 

The location of individual frequencies is governed by 
the asymptotic equation (Tassoul 1980 ): 



Av n+ -Z + e 



5u { 



'01- 



(1) 



In this equation, v n j is the frequency with radial order n 
and angular degree I (neglecting rotational splitting so there 
is no dependence on m), and Av and Svoi the large and 
small separations while e is an offset coming from second 
order and near-surface effects. This theoretical relation is 
used to identify the individual modes, from the smoothed 
power spectrum. The entire procedure is illustrated in Fig- 
ure Q] 

We have determined 14 (I = 0), 13 (I = 1) and 16 (I = 
2) frequencies, which constitute a total of ten more frequen- 



cies than the initial asteroseismic analysis ([Christenseri- 



Dalsgaar d et al.p OlO). The frequencies were determined 
with an average standard error estimated to be cr v =0.75 
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Fig. 2 An echelle diagram showing the observed frequen- 
cies (red) and the best model (blue). 



/iHz (compared to a previous accuracy of 1.4 /iHz). They 
are illustrated in Figure [2] 

2.2 Modelling 

To calculate evolution models and adiabatic oscillation fre- 
quencies, we use the Aarhus codes ( Christensen-Dalsgaard 
2008a b). Evolution tracks were calculated for a grid of 
models with mass range 1.33-1.70 M and metallicity 
0.015-0.035. The convective core was modelled with over- 
shoot in units of the pressure scale height of aov = 0, 0.1 
and 0.2. Based on the observed large separation an estimate 
of the density is obtained ( W hite et al.|2011 ), which is used 
to select relevant time steps in the stellar evolution tracks, 
for which frequencies are then calculated. 

The calculated frequencies are compared with the obser- 
vations, after correcting for near- surface effects that are not 
modelled properly. The empirical correction law is given by 
( Kjeldsenetal.|2008) : 

^obsH l6 



^obsO) - ^best(^) = a 



"0 



(2) 



where is a reference frequency and a and b are parameters 
to be determined. In addition, we set ^best(^) = ^ref(^)> 
where r should be close to one as the reference model 
should be close to the initial best model. In practice, we 
set b = 4.90 (the solar value) and we find the best fit for r 
and a (for each model individually). We can then calculate 
the reference model, for comparison with the observed fre- 
quencies. The best model is then picked by comparing the 
observed frequencies with the corrected model frequencies. 
We use a y 2 minimisation to determine which model fits 
best: 

1 „ / ^ ,„^2 
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The frequencies of the best model are shown together 
with the observational frequencies in the echelle diagram 
in Figure [2] With the exception of a number of frequencies 
(between 1200 and 1400 /iHz) for I = 1, the frequencies are 
modelled very well, resulting in x 2 = 3.65. Gaussian noise 
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Table 1 Parameters obtained for star (HAT-P-7) and its 
planet (HAT-P-7b). 

HAT-P-7 



Mass (M ) 


1.361 ±0.021 


Z 


0.01904 ± 0.0015 


[Fe/H] 


0.13 


Age (Gyr) 


2.19 ±0.12 


Radius (R ) 


1.904 ±0.010 


Temperature (K) 


6259 ± 32 


Density (g/cm -3 ) 


0.2781 ± 0.0017 


Luminosity (L ) 


4.996 ± 0.098 


\ogg (cgs) 


4.98 ±0.13 


Ellips. var. d (ppm) 


59 ± 1 


HAT-P-7b 


Orbital period (d) 


2.20473506(11) 


Orbital inclination i (°) 


86.68 ±0.14 


M p (Mj) 


1.741 ± 0.028 


Rp (Rj) 


1.431 ±0.011 


Tday (K) 


2470 ± 10 


Tonight (K) 


1772 ± 10 


A* 


0.193 ±0.002 



was added to the observational frequencies to create 1000 
different frequency sets. For each, the best model is calcu- 
lated and the stellar parameters are calculated. The average 
and standard deviation of these parameters are presented in 
Tableffl 



the bottom of the transit depends on the limb darkening, 
modelled by a quadratic intensity profile ( |Winn|20T0 ): 



I(X,Y) oc 1 

vT 



ui(l - aO - ^2(1 - m) 2 



(4) 



with fi = y/l — X 2 — Y 2 and (X, Y) the coordinates on 
a unit circle. The coefficients Uk are constants defining the 
precise shape of the limb-darkening law. We have simply 
used them as fitting parameters. The modelled transit is 
shown in Figure [3] We have also made a new estimate of 
the planetary mass, using the radial velocity values derived 
by |Pal et al.|p008) . 



3.3 Phase variations 

The planetary system is tidally locked, which causes the star 
to always 'see' the same side of the planet, while a distant 
observer sees the planetary dayside and nightside at differ- 
ent times. However, for HAT-P-7 the situation is slightly 
complicated. The close orbit of the planet has induced a 
tidal distortion on the star: rather than being oblate due to 
rotation, it has its longest axis towards the planet and the 
shortest axis perpendicular to the orbital plane. This causes 
ellipsoidal variations (which peak at phases near 0.25 and 
0.75), which should be added to the phase variations ([Pfahl 
eFaLl[2008l |Welsh et al.pJIO) . HAT-P-7 is the only cur- 
rently known system with ellipsoidal variations caused by a 
planet. We model the flux throughout one orbital phase as 



3 Planetary properties 

3.1 Period 

The time of each individual transit is measured and used to 
determine the planetary period. Assuming a perfectly Ke- 
plerian orbit, we fit a straight line through all transit times 
(taking gaps into account; whenever a transit is measured 
incompletely it is ignored). Due to the large amount of mea- 
sured transits (286), the period is determined very accu- 
rately. An O-C plot has been calculated, which compares 
the observed and the calculated transit times to look for a 
gravitational influence from additional bodies in the sys- 
tem. No sinuisoidal pattern was detected, and variations are 
at the level of the observational error level for individual 
transit times (« 0.00015 days). After the period has been 
calculated, the time series is folded (calculated modulo the 
period). 

3.2 Transit shape 

From the folded time series, the transit is modelled. The 
transit duration determines the impact parameter, which is 
directly related to the orbital inclination. The transit depth 
determines the relative planetary radius, which can be com- 
bined with the stellar radius (derived in the previous sec- 
tion), to obtain the absolute planetary radius. The shape of 



/(*,<) 
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cos(2a) 

7T 



sin(a) ± (tt — a) cos(a) 



(5) 



with the orbital phase & G [0,1] taken as at maximum 
radial velocity of the star, so that the phase angle a is de- 
fined as cos a = — sin i sin27r<£. The first part of f(<&,i) 
is the Lambert law for a sphere (Charbonn eau et al.|1999| ), 
where A g is the planetary albedo and a is the semi-major 
axis. A non-zero value for d corresponds to the ellipsoidal 
variations. We refer to Table[T]for values and to Figure[3]for 
the best fit. 



3.4 Occultation 

When the planet moves behind the star, an occultation oc- 
curs. The depth £fl ux (A) gives information about the plan- 
etary dayside temperature. By combining the information 
with the phase variations, we can also obtain the night- 
side temperature. We approximate both star and planet as 
a blackbody and use 



^flux(A) = 



(6) 



Exp*)' 

where B\(T) is the Planck function (Winn 2010), which 
we integrate over the Kepler bandpass. The exact occulta- 
tion depth is calculated by simply comparing the averages 
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Fig. 3 After the planetary period is determined, the obser- 
vations (blue) are folded. This allows for a detailed study of 
the planetary transit (top figure), phase variations (caused 
by planetary day-nightside and ellipsoidal variations, mid- 
dle figure), and planetary occultation (bottom). The black 
line shows a data fit, the colored data points are used in cal- 
culating the occultation depth. 



of data points inside and just outside the occultation, as 
shown in Figure [3] We find an occultation depth of 71.85 
± 0.23 parts per million. We assume that all planetary flux 
comes from radiation (albedo zero), so the quoted temper- 
ature value is an upper limit. The indicated accuracy of the 
quoted temperatures should be seen as the internal error of 
using a black-body spectrum to describe the spectrum of the 
exoplanet. Since the black-body flux in the Kepler bandpass 
for exoplanet temperatures will change a lot even for a small 
change in temperature one can estimate the black-body tem- 
perature precisely. We find the occultation to be 4.90 ± 0.25 
ppm deeper than the value just outside transit. 

4 Discussion 

The asteroseismic analysis of the host star has followed 
largely the same methods as applied by Christe nsen-| 
Dalsgaard et"aT] ( |20T0| ) for Kepler's Q0 and Ql data. With 
eleven quarters now available, all parameters are determined 
to higher accuracy. The agreement is fairly good, but we 
have derived a slightly lower mass, radius and metallicity. 



We have not been able to model a few frequencies prop- 
erly (see Figure [2]), a problem that did not present itself in 
the earlier analysis by Christensen-Dalsgaard e t al. ( 2010 ) 
because of the fewer frequencies observed, and a lower ob- 
servational accuracy. Ellipsoidal variations were found by 
Wel sh et al. ( 2010| ), and we confirm that the phase curve 



cannot be fit without including this effect. | Welsh et al. 
( |2010| ) find a variation of 37 ppm (no errors quoted, based 
on Kepler Q0-Q1 data), which is considerably lower than 
our value of 59 ± 1 ppm. All planetary values have im- 
proved accuracy. Most are in general agreement with earlier 
work. Our planetary temperatures are significantly lower 



than all earlier temperature determinations ( [Borucki et al. 
[20091 IChristiansen et al.|20T0l|Pal et al.|2008U Welsh et al 



2010). They all used an occultation depth of 130 ± 11 ppm 
( [Borucki et al.|2009] ), which is almost twice as high as the 
value derived in this work; the analysis was done based only 
on Kepler Q0 data, and the modelling did not include ellip- 
soidal variations as they were not yet discovered in HAT- 
P-7 at the time. On the other hand, the temperature in this 
work assumes a blackbody, which is an oversimplification 
and also causes the temperature error to be underestimated. 
|Welsh etal.1 ( |2010| ) and |Christiansen et al.| ( |2010| ) used more 
realistic atmosphere models, and it should be interesting to 
see how they are affected by the revised occultation depth. 

This work shows the excellent quality Kepler can now 
deliver, both for studying stellar pulsations and planetary 
transits. Asteroseismology is likely to take up an increas- 
ingly important role in studying host stars of exoplanets. 
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